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TWO APPROACHES TO GRAVITY REFERENCED ORIENTATION 
aa GRLENTATION 


dey First Approach 


Given camera parameters and overlapping exposures of an object the 
shape and configuration of an object may be determined by the use of 
conjugate images in suitable reciprocal collinearity equations, If only 
the shape of the object is required all exposures may be referred to any 
orthogonal coordinate System inherent to the object, However, if the 
object's orientation is a function of the direction of gravity then an arbitrary 
coordinate system with an unknown relation to the direction of gravity 


is unacceptable, 


In the absence of any gravity referenced control data the search 
for the gravity direction in the arbitrary coordinate System must be 
found in those motions under the influence of gravity. The motions of 
the camera associated with pitch, roll, and yaw are a coupling consequence 
of the drag resistance of the camera through a space under continuous 
restoring force of gravity. Translations from camera station to camera 
Station are mechanically imposed while the changes in orientation occur~ 
ring during translation have the direction of gravity as the axis of 
equilibrium. These axes of equilibrium are illustrated in Figure 1, If 
Xe Vy and Zz, are the axes of a camera vehicle S is the direction of 
gravity and Y_ is the horizontal direction of motion. At rest x» Ye 
z,, coincides with Xo vo and Zee Pitch (t) is the vertical angle ve 
axis defines with the horizontal as a plane, Rent (9) is the rotation 
of Xe 2 plane about the ae axis from the = ag plane. Thus both pitch 
and roll refer to the ne te plane but 7 is a vertical angle that always 
contains a and p is not a vertical angle except when 7 = 0, p lies in 
the x» 4 plane when 1 #0. Yaw (c) is the rotation about Z. of the vi 


axis or the horizontal angle the vy, axis makes with ue axis, 


Assuming the camera to be rigidly attached to its transport vehicle 


the camera axes may be treated as vehicle axes with constant angles relating 
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7? *POLL 


IN 


WT = PITCH 


Figure 1 PITCH, ROLL AND YAW 
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the vehicle axis to the camera axes, Let To» Pos and BG be the constant 
pitch, roll, and yaw of the camera with the vehicle at equilibrium whence 
changes in vehicle orientation do not alter the camera relation to the 


vehicle. 


It is believed that a given vehicle had undergone normal pitch and 
yaw rotations with minimal or nominal zero roll rotations. This means 
the x axis lies in the ae ty plane and defines a level line which with 
horizontal change associated with yaw defines a horizontal plane whose 
perpendicular is the direction of gravity. This geometry is shown in 
Figure 2 with two different orientations. It is evident if yaw is a 
constant x does not change its horizontal direction without which the 
direction of - cannot be determined, Figure 3 illustrates the proposed 
solution based on the assumption of significant pitch and yaw and small or 


negleetable roll. 


Assume p = O and ee or the axis of pitch is a level line. X' 

in the presence of yaw about Z, defines -a plane whose perpendicular is 

eS X> Y Zo is the arbitrary static object space coordinate system to 
which the camera undergoing yaw and pitch rotations refer. The camera 
orientations with respect to the arbitrary datum is established. It is 
desired to refer the camera orientation to the direction of gravity and 

aa ae system normal to A ae and ae lie in a horizontal plane but 

because of yaw are different for each exposure. However, the camera axes 
defines constant angles x x", y. X', and Zo X' with the yawing a (x) 
axis. It is necessary to determine both these constant angles and the 
instantaneous direction of x" in order to determine the direction Zo 
Given the orientation matrix of a camera axes referred to the arbitrary 


datum, 
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‘Figure 2 CHANGE IN 7 AND JY WITH Y =0 FOR TWO EXPOSURES. 
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HORIZON PLANE Y’' 


Figure 3 GEOMETRY OF SOLUTION FOR DIRECTION OF 
Zg IN Xa,Ya,Za COORDINATE SYSTEM, 
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As usual the columns are the direction cosines of the camera x, y, z axes 
with respect to the Xo Yo Zo axes and the rows are direction cosines 

of the Xo Yo Zs axes with respect to the camera x, y, z2 axes. If one 
knew the direction cosines of the instantaneous a axis. referred to the 
arbitrary datum one could write for the constant camera-vehicle angles. 


the following equations: 


COS Oy, COS A, + cos By cos BS. + cos ¥yr cos y, = cos xox! 
COS Oy, COS a + cos Byr cos By + cos Vyu cos Yy = cos yox' 
cos zoX' 
+ + 
COS Oy, COS O, cos Byr cos B COS Yy1 COS Y, 


Since the right hand terms are constant 


cos xoX' cos a, + cos yoXx' cos oy + cos #oX' cos a, ™ cos o 
cos xoX! cos B, + cos yoX' cos B. + cos zoX' cos 8. = cos Byt 


cos xoX' cos Y, + cos yoX' cos ie + cos zoX' cos Y, = C08 Yyr 


Now since ay is perpendicular to ae 


cos oO z. + cos Bye cos Bo + COS Yy+ COS Y, = 0. 


8 & 


While Gyre Byrs and Yer Vary for each exposure they are expressible 
in terms of the known arbitrary orientation matrix and the constant angles 


of the vehicle level X' (x) axis with respect to the camera x, y, z axes. 


Substituting 


(cos a cos xoX') cos a. + (cos ae cos yoX') cos oy + (cos G, cos zoX') cos 


& & & 


Z 


Qa 


“2 


+(cos By cos xoX') cos B. + (cos By cos yoX') cos Pe + (cos 8, cos zoX') cos 8. 


& 8g 8 


' ' 3 1 
+(cos Yz cos xoX') cos Ye, F (cos Yz cos yox ) cos % + (cos Vz coszzoxX )} cos Ne 


& g 8 
= QO, 
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Division by cos Yz, cos zoX' gives 
g 


(tan "2, tan Ny) cos a, + (tan "2, tan Ex) cos - + (tan "2 COs o, 


+(tan 2, tan Nyt) cos B. + (tan ae tan Eyr) cos By + (tan ae cos 8 


+(tan Nyt) cos y, + (tan byt) cos ie = - cos Me 


which contains 8 unknowns in parentheses, Assume an 8 x 8 is solved 


to obtain preliminary values of tan Ny » tan by » tan yrs and tan Ey 


g§ 
Then 
tan Ny 
cos a, = g 
£ (1 + tan Ny + tan by ) 1/2 
& g 
tan 7 
cos Bo = g 
g (1 + tan ny + tan by ) 1/2 
g & 


cos y, = 1 
g d+ tan7n, + tan7é, be 


& & 


tan ny, 


(1 + tann,, + tan7E.,) 1/2 


eos xoX! = 


tan byt 


cos yoX' = 
(1 + tan7n,, + tan7é 1) 1/2 


cos zoX' = 


1 
(1 + tan7n,, + tan7é, 4) 1/2 


7 
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Inasmuch as op = 0 
flight will yield poor 
get good results 40 to 
by Least Square values 
the right hand term of 
residuals obtained are 


90° or cos (90° + —) = 


any 8 x 8 formed from sequential exposures in 
results because of the small changes in Eo To 
50 exposures must be normalized to an 8 x 8 where 
may be obtained and residuals generated. Since 
the original equation is zero or cos 90° the 

the cosines of angles less than or greater than 


+ sino. Thus the residuals are a consequence 


of the assumption 9 = 0. A large number of condition equations will 


include the necessary spread in yaw - especially if the values are 


selected from different flight directions. 


Assuming those values of sin p having 3 times the mean deviation 


from the average sin p 


are rejected from the solution an iterative form 


of the equations is normalized and solved. The iterative form is a 


4x4 


A Any +B mee C 


where 

A = (tan n'yr cos 

B = (tan niys cos 

Cc = (tan n'y cos 
& 

D = (tan n'y cos 
& 


The primed yn" and 


normalized 8 x 8 which 


Any + D Abyy = sino 


t + 2.,% 
o + tan & x1 cos a, cos a) sec’n', 
+ : + 2e% 
B. tan & x1 608 By cos B) secté 
+ : + 2n' 
O tan & z cos By. cos ¥5) secn' ys 


' 2rt 
oy + tan & z 608 e + cos yy) sec“é x! 


&' denote the best values obtained from the 


are still considered preliminary. Solving 


normalized equations for Ano» MEye At» and Agyy 


1 
n z + tan, 


& & 
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= ct 

fy TE gy F BAY 
& & & 

Ay > n'y8 + TANS 

ey Et ye + URE ys 


Naturally values tan n'y » tan &' 


g 
with each iteration until the © sin* p = a minimum, 


y » tan Nyt and tan nyt are revised 
8 


An important point not mentioned is the fact that the cameras cannot be 
interchanged nor can the camera x, y axes be different from exposure to 
exposure, Each camera defines different constant angles with the vehicle 
and any one camera x, y axes defines constant angles if the camera x, y 


axes are not rotated from frame to frame. 


Suppose for some reason the camera x, y axes were rotated 90°, 180°, 


or 270° as is indicated in Figure 4. 


In each case for convenience in relative orientation the camera 
axes has been rotated from the original orientation Xa» Y, to x, y for 
which the matrix is given. It is only necessary to restore the matrix 
to the original axes when the rotation is 90°, 180°, or 270° to make the 


sign and/or column change indicated below: 


90° Rotation 


a a a —e. oe ae 
@) Yo 20 y 

by. OB b _|-> Ba, RY 
oO %o 26 y 

ee c Ss as ce ae 
[@) Vo Oo y 


9 
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90° Rotation 


Xo 


180° Rotation 


270°Rotation 


Figure 4 ORTHOGONAL ROTATIONS 
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180° Rotation 


a2 a a) = -a a 

oO %o ie] y 

bb Bes eter. be. “Be 
fe) Yo fe) y 

Cx °. co ~ey ia Pa 
ie) Yo O° y 


a a a a 7 ee ao 

ie) Yo 25 y 

ee en ee 
1°] %o °o y 

Cy Cc c. c ie e. 
[e) Vo [e) » 


A test solution of an 8 x 8 for each of two cameras yielded both 
éncouraging and discouraging results. The standard coordinates tan Ny 
and tan Eo were roughly the same in magnitude and sign for both cameras. 
This would not be possible unless the solution was based on correct 
theorems. The values tan n'y and tan byt showed large variations as might 
be expected from 8 closely grouped exposures and/or containing an unknown 
camera x,-y rotation. The signs of tan Tyr and tan byt were opposite 


for the two cameras as was expected. 


The consistent values of ng and by and the consistent signs of the 
Nyt and byt values suggest the solution is valid and the assumptions are 
valid, The large spread in the values of Nyt and byt obtained by two 
camera solutions indicates the x, y axes must be rotated back to the 
original camera orientation and a large over determination is necessary 


to insure a spread in yaw and to average out the variation of 9 from zero. 
In view of the sensitivities of the 8 x 8 and the linearized 4 x 4 


to both unknown x, y rotations and departures from the assumption of p = 0 


an unduly large effort is necessary to pre-process each orientation matrix 
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employed. To surmount this difficulty a solution is proposed not dependent 


on 0 = 0 and independent of x, y rotations. 


2. Second Approach 


It may be noted that the z columns of the orientation matrix are 
unaffected by an x, y rotation. The following approach is suggested to 
provide a solution that is not dependent on p = 0 and is independent 
of an x, y rotation. The 8 x 8 and linearized 4 x 4 indicate p is not 
guite as small as thought and 7 is not quite as large as thought. If 
this is true, each axis (x, y, z) of each camera with nN exposures 
generates a right cone with a sinusoidal surface and whose axis is the 
direction of gravity. An equation is written for each exposure of each 


cameras 


x axis 


eon xe cos a 4-{ C08 Bo cos 6 +f O8 Y¢ cos Y 
cos x 1 cos x 1 cos x 1 
& & g 


ce cos a +/ 508 Bo cos B +{ °° Yo cos y = 
So x Saeed x ee x = 1 
cos Xp 2 cos K, 2 cos X 2 


fl 


8 


eee eee mm eee ees e eee ee ee es Eee eeeHeEDn eer eee ORE EeEH eR Odo eo eee een dnEene 


Fs ie cos a {= Bo ] cos B SS ife } cos y =1 
x a x sis 
n n n 


cos Xy COs X, COs Ey 
y axis 
al os Bo ) cos B Ye cos y =] 
cos Ve Yq cos oe Yq cosy, Mat 
= 1 


Hoe NG cos ao. +f °°8 Bo cos 8 + re a cos 
2 \cos ie 2 \cos Te V9 


ee eeeeeeaereeeeceereeerneeeesereeeee Heese eeeereereseenereeeeerereee 
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dao € cosa + °° Bo cos B +{°°8 Yo cos y = 1 
cos vp Yn \cos Ye Yn \cos y Yn 


& 
zZ axis 
sedeas ¢ cos a +{ 608 Bo cos 8B Os ean (6! cos y = 1] 
COs Z 1 cos Z 1 cos Z ey, 
g g 
hie we cos oF ae Bo cos 8B (CoP Ve cos yY = 1 
n Zz Zz 
cos Zz cos Z 2 cos Z 2 
& g & 


Coote oreteoeeseervneeereeeeeenerneaeseeenenaeeeeeeeseeeHesee Hes ePeeaeaeenenenee 


eee “Gg cos &, EN aeies Bo cos Bo 4( C08 va cosy, 2 y 
cos 2. n cos Z cos 2 n 


& 


The geometry of the above three solutions is illustrated in Figure 


g" “Ss 
direction and gravity. Cas Bas and Yg are the direction angles of the 


5. X y y. and Ae are the cone angles each axis generates with the 
direction of gravity in the arbitrary coordinate system. 


Only the third array is independent of an x, y rotation. However, 
if there were no x, y rotation the first and second array would yield 
the same direction of gravity except for the cone angle generated, The 


approximate orientation of the system at rest is 


= 114° (65° 
x, (65°) 
= 0° 
Yg 
Zz = 24° 
g 


This means that the x solution in the absence of an x, y rotation would 
yield results comparable to the z solution whereas the y solution being very 
approximately a level line (direction of motion) would yield a poor 
solution. The y cone degenerates to a sinusoidal plane. Since the 
departures from cosine 90° are equal to or greater than several degrees, 


the cosine of 90° + 2° has twice the numerical magnitude -as the cosine 
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Figure 5 SINUSOIDAL CONES GENERATED BY x,y and z AXES ABOUT Zg 
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of 90° + 1° for example. cos Vy being a denominator can affect the 
numerator cos Oa cos Bae cos Ye by several 100% with an error of several 
degrees. Assume the values in parentheses have been solved from a 


normalized 3 x 3 


cos a = 1 
Se PR TN ee PURE, Ne: 
COS Og + cos Ba +{ 68 Ye 1/2 
cos xX cos x cos xX 
8 g g 
cos Ve 7 1 
$$ $$$ 
le a + (= “| ‘| a8 NG | 1/2 
cos cos cos 
Vg Yg %g 
cos Ze = 1 
a eS eee Ae od Oak 
cos GO, 4 cos Ba 4{ COS 'Y¢: 1/2 
cos Z COs Z cos Z 
8 g & 
whence 
COS OU, = us ne cos xX_ = ces G cos y_ = eG cos Z 
cos X cos y & \cos 2 & 
8 g g 
cos 6, = ae G cos x_ * a G cos =| 208 Be cos Z 
GO [Fs 8 Yg aa a g 
cos xX cos y cos 2 
cos Y, = ROR TG cos xX = caueeeuce cos y_ = eee te cos Z 
Os re SY. See Bo eeaoe. ss g 
cos x cos y cos Z 
8 g 8 


The particular values of ap y and 25 for each exposure are obtained 


gs? 


by substitution of the equated values of cos o cos Bas COS Yq back 


into the original condition equations: 


cos a, cos a + cos cos + cos cos = cos x 
G x Bg B. Yg Vx n 
n n nl & 


15 
Approved For Release 2001/08/02 : CIA-RDP80T00703A000800040001-0 


Approved For Release 2001/08/02 : CIA-RDP80T00703A000800040001-0 


cos Ga, cos a + cos 8, cos 8 + cos y, cos y =;cos y 
G vy G ys G Ye a 
cos vel cos O, -- cos Bo cos Bo + cos Yo cos Y, = cos Zo 
n n n. g 
whence 
(pitch) An, a. 29 Te 
g oO 
(tilt) At = zo % 
n g By 


(roll) ap = Py ~ Py 


where 
sin x 
tan op = Bo 
0 eareienree ra 
cos Z 
& 


tan o 


These quantities are only of academic interest, 


To test the above equations, the x, y, and z columns of eight 
camera orientation matrices were normalized and solved by the method 


of Least Squares. The results are given below: 


“x y Z 
COS Ol, 0250543 -.1388537 - 1172539 
cos Bo - 2918100 - 9568414 - 2801572 
cos Yq -9561481 -. 2552922 -9527663 


As anticipated, the direction cosines of x and z are consistent 


to the first order where as those of y.are absurd. By treating cos ve 
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as zero the solution for cos as Cos Bos and cos fe is as follow: 


cos a, = -0962929 
cos Ba = .4908134 
cos Yq = 8503334 


which means treating y as sweeping out a sinusoidal plane in place of 
a conic removes the absurd aspects of the values even though the results 


are not as logical as those obtained with z, 


Were it not for the problem of x, y rotation a solution embracing 


a conic could be employed with the rows as well as the columns. 


X row 
ona cos + rage cos Sates cos a =1 
eects (COR ety ee oo) pecs -seee : 
cos a cos oO cos 4 
G G G 
or 
Y row 
cos x cos y COS Z 
g cos + g cos 7 + g cos = 1 
cos 8G cos BG cos Bs 
or 
Z row 


il 
| aed 


cos x cos y cos Z 
g cos z + g cos + g cos 
cos Ye y cos Ye 


Solution of normalized 3 x 3 in each of the above forms yleided 
reasonable results for the third row but for the reason of weak geometry, 


absurd results in the first and second row, 


The evidence points toward the most reliable solution being obtained 


with the third column. To this end 67 starboard and 71 port third column 


LF 
Approved For Release 2001/08/02 : CIA-RDP80T00703A000800040001-0 


Approved For Release 2001/08/02 : CIA-RDP80T00703A000800040001-0 


matrices were employed in a Least Square determination of cos Qos Cos Bas 


cos Yq and cos eae 


Starboard (67) 


Cosine 


Zoe - 9072857 


oF » 1948335 
Be - 2858874 
Ye 9382475 


Angle 
24° 52! 
78° 46' 
73° 23! 
20° 14' 


The results are given below: 


Port (71) 
Cosine Angle 
» 9085228 24° 43! 
- 1669294 80° 23! 
» 3628651 68° 43! 
-9167680 23° 33! 


Inasmuch as the interlocking angle between the two cameras is 


some value greater than 48° and inasmuch as: the bisector of this angle 


approximates the vertical, solved for values of z 


Starboard 24° 52' and 


Zz, port of 24° 43' establishes confidence in the above direction angles 


of the vertical. 


The difference in the starboard and port gravity angles 


may be attributed entirely to the noise of the sinusoidal assumptions. 


3. Precision of the Gravity Direction Angles 


In order to demonstrate the precision of the determination a 


numerical example of 8 starboard 3rd column values normalized to a 3.x 3 


is given: 


(= =| (= Bo 

oe cos + | 

cos Z Zz cos Zz 
& g 


a 


«2978860 
-.2073995 
-.2193539 

- 1588703 
-.0026143 

- 1853262 
-.2711262 

4051708 


+9787331 
«4288371 
+ 4242990 
- 6334470 
0984186 
«5969892 
«3807350 
0442197 
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COS Yq 


cos + (sate 
Zz cos Zz 
& 


Cc 


«7591650 
«8792521 
+ 8785523 
«7573012 
9951417 
- 7805498 
» 8840427 
»9131710 


el oo = = 
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Normal Equations 


abe © [aa] + ape Bo [ab] + oon kG [ac] 
cos Z cos Z cos Z. 


[a] 


& & 
( eer 7 [ab] i cos: Bi [bb] i ron Gg [be] = [b] 
cos Zz cos 2 cos Z 
& g g 
picts © [ac] ‘t pee BG [be] S ive oe = [c] 
cos Z cos Z cos Zz 
g 8 g& 
The brackets denote sums. 
24771323 . 1166048 »2437452 = 0.3467600 
1166048 1.6131149 2.4138944 = 2.9888415 
» 2437452 2.4138944 5.9097528 = 6.8471758 
Solving by Crammer's Rule 
(= we \ . 1267239 (= fe = 1.0297208 
cos Zz cos Z 
& 8 
i BG 3027854 
cos Z 
g& 


cos 25 = lL 
Tae en Pane RNG. Pame. ae ee 
| “s) + (= =) y (= te| i= 
cos Z| cos Z cos Z 
& & & 


= 0.9252666 
$e 225 175 


= = ° ' 

COS O, = 1172539 oval 83° 17 

= = ° ' 

cos Ba = .2801572 Ba 73° 44 

= = ° ' 

cos Y, 9527663 Ye 17° 41 
19 
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To determine the residuals the values of cos dg » cos Bas and cos Ye 


are substituted in the original equations: 


cos ze - (cos "a cos Oo + cos an cos Bo + cos ‘2, cos Yq) 


= cos 2 —- COS Z 
8 n 


where Z. is the Least Square value and es is the value obtained from any 


condition equation n ° 
now cos Z = cos z cos A. + sin z sin A 
Gy & a & Ze 
a a 
To the first order cos A = 1 therefore cos z =- cos Zz =- sinz.A 
Z & g g 2 
. Bn n 8 
cos Z —- cos 2Z 
g g. ; 
he = 
PO a oa z! (in minutes of arc) 


sin z - sin l' 
g n 


With this formulation the following residuals were obtained: 


Bo is “HPO 
81 

A, mw 75.0 
Bo 

An. = 44.8 
83 

An =) (69,3 
By 

An =) (45.4 
24 
B 5 

Be. tee OTL 
86 

An =) 73.5 
o 
37 

42, = -42.3 
3 g 
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The values are tilt oscillations on the surface of the cone 
generated by the starboard optical axis and are equal in magnitude to 
what might be expected, Since the tA = 3'.5 the value of z obtained 

g 


is the most probable value exhibited by the 8 starboard exposures, 


Thus 


= 57'.7 (ignoring signs) 


where n = the number of equations and u = the number of direction angles. 


In order to determine the mean error of Cas Bas and Ye it is 
necessary to compute the relative weights of the gravity direction angles 


from a modified form of the normalized 3 x 3. 


COS OG [aa] + cos Bo [ab] + cos V8 [ac] = 1, 0, 0 
COS OG [ab] + cos Bo [bb] + cos Yo [bc] = 0, 1, 0 
COS O, [ac] + cos Bo [bc] + cos Yo [ec] = 0, 0, 1 


Using the first right hand column for Pas the second for re and the third 


for NG the following equations and numerical values are obtained 


Po, = A = _1.7293830 = .4666159 
[bb] [cc] — [bc]? 3.7062241 
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P B = A = 1,7293830 = .6265149 
[cc]-[aa] - Lac]l* 2.7603222 

ees A = 1.7293830 = 2,2873243 
[aal.[bb] - [ab]2 ~.7560725 

where 


A = [aa {[bb]-{cc] - [be]?}- fab] {[ab]-{cc] - [be]-[ac]}+ lac] {fab]-[be] - [bb] -Lacl} 


An alternate method of computing weights was employed as a check. To 
determine the weight of variable a? 8a and Yq are expressed in terms of oS in 
the Be and Ye normal equations. These values are substituted in the normal 


equation a The resulting coefficient of Qe is the weight of a To determine 


G° 
the weight of variable 8 


G* 
q? % and Yq are expressed in terms of Be in the Ob 
and Ye normal equations. These values substituted in the normal equation in 


Ba form the weight. of Be in the resulting coefficient. Similarly the coefficient 


of Yo gives the weight of Yo when o and B. expressed in terms of Yq are 


& 
substituted in the normal equation in Yer 


P o 
g 
[bb] cos G + [be] cos a -[ab] cos de 
[bc] cos G + [ce] cos er ~[ac] cos Q& 
cos Ba = 


(test [be] - [ab]. Lac]\ ogg g 
[bb]. [cc] —- [bc] G 


COS Yq = [bc]. fab] - [bb]. Lacl\ cos a 
[bb]- [ec] = [bc] c 


Substituting in the qa, normal equation 


G 


P Oo= [ab] + oat ac]. [bc] - [ab]- [cc , + [ac] (ES - pe lee 
[bb]. Lec] - [bc] [bb]-[cec] = [bc] 
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P Ba 
[aa] cos Ota + [ac] cos a -[ab] cos Ba 
[ac] cos ao + [cc] cos G = -[bce] cos Be 


cos a 


i?) 


*{ [be]-[ac] - [sb] sLeel } eos Bo 


[aa]-[ec] - [ac] 


<= : [ac]-[ab] ~ ibe] -Laal | cos B 
eG { [aa]-{cc] - [ac] 


P g.= [ab] [bc]: [ae] - [ab]*[ec]|+ [bb] + [be] { [ac]+[ab] - [be]: [aa] 
: faal’ [cc] - tac] Taa}-[ec] - Lac] 


PY¥¢ 


[aa] cos Oo + [ab] cos Bo = -[ac] cos Y¢ 


[ab] cos Oe + [bb] cos Be = -[bc] cos Ye 


COS Og = [bce].[ab] - tac}. [bbl | COS Yq 
[aa]. [bb] Pa [ ab | 


cos 6, ={_[ab].[ac] - [2a] {bel } cos Yo 
G Taal. [bb] — [ab] 


= be}.{ab] - [ac]-[bb] | + [bc]} [ab »{ac] - [aa]-[bc] [| + [cc] 
PG eee {Pe}-fet = taste _ fh te aalefbb] —~ [abl- 


Exactly the same numerical values were obtained P al 4666159 
13 Ba = .6265150 
P Yq = 2.2873241 
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The mean error of each direction angle is computed with 


= = = ' gn ° 
ed, = _-g 15.0 109'.8 = 1° 49,8 
YP de ‘ 
@ 
e 8 = 79 = io = 94',8 = 1° 34,8 
YP 8 ; 
fe 
ene 
— = = ’ 
2 ¢ 4 8 , 75.0 49'.6 
VP Ye 7 


While these errors by ordinary standards are large those based 67 on star- 
board exposure and 71 port exposures are much smaller. Even these error values 
are adequate for the purpose of referring the object space coordinate 


system to the direction of gravity. 


4. Further Adjustment 


Inasmuch as cos Go» Cos Bas and cos Yq are dependent functions 
it would seem that improved values would be obtained if the condition 
equation were linearized in terms of two independent Eulerian angles. 


Let 


cos O. = cos 9G sin Yo 


cos Bo = sin 8G sin fe 


COS Yq = Cos Ye 
whence 


(cos 6, sin YQ) cos a, + (sin 8G sin Yq cos B + cos Yq COS Y, = cos a 
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If 


7 fe) + cos f s + Cc = -si 
GOs) @ (cos a, cos a, + cos B, cos 8B, + cos y, Cos see) ae AL 


& 


G 


or 


A.AO,+B, AY@ =e of 


G Zz 


g 
where 


A = (cos By COS O, + cos a, cos Bb) / sin Ze 


G 


to 
i) 


[ (cos tg cos a, + cos B, cos B) cot y, - (cos Yq cos ‘2 tan YQ] / sin a 


n condition equations of the above form may be normalized to a 2 x 2 in AGG 


and AYqs A test solution demonstrated that the best values had already been 


obtained inasmuch as on the first iteration Ad, and AY were significantly 


G 
zero. Had the condition equation been linearized, it would have the following 


form: 


sin O, COs a, Ada + sin Bo cos B. ABG + sin Ye oS Y, Ay, = sin ar A. 


G 
& 


Assume there is a particular value of as Bas and Ye that satisfies the equation: 


+ + - Z =Q 
cos el cos a, cos BG cos BC cos Ye cos wy cos g 
n n n n n n 
Since cos Z —- cos Zz = -sin Zz A 
& Bn & 20 
n 
= + Ss 
cos oe cos Oo A co Qo 
n n 
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cos Bo = cos Ba + A cos Ba 
n n 


COS Yq = COS Y, + Acos Ye 
n n 


and 


A cos BG = cos Bo - cos Ba = -sin Ba AB 


A cos Ye = COS Yq — cos Ta: = -~sin Ye AYa 
n n 
Substituting we obtain to the first order the form obtained by linearization: 


G 


+ As + 
A cos i cos a, Arcos Pa cos B. A cos Ye. cos y, 


= cos Se - (cos Oe cos a, + cos Bo cos Be + cos Yo cos Y,) = -—sin Ze a 
or + 81 4 cos a (Aa.) + sin 86 cos B_ (AB) + si0 Y¢ cos y, (Ay,) = +A 
: a G rs Zz G aac! zZ G Z 
sin z sin 25 sin 2. g 


New normal equations in Ada» ABa » and AY@ are easily formed since all 
coefficient changes are constant. It is necessary to accept the values of 
Qos Bas and Yg as the best possible if the correction determined from the 


revised normal equations are significantly zero. 


sin2a sin %. sin B [ab] AB 4 sin a, sin y, Padi awce Packs 9 sin 
aoe aaa ae mat G fap) ORs sin* z aN ae ears 
g sin Zo g & 2 
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£ Sarge 


B 


sin 4, sin 8 2 i | si 
ae — [ab] do, + 818 8o fbb] ag, + SD SG AIR NG [be] aye =[b.a, 1 ay 
aoe 45 sin’ z sin Zz g sin 2 
g g & 
2 : 
: : in Y sin 
sin % 810 Yq [ac] ho, + sin 8,.810.¥, [pe] a eG [cc] Ay,= [cA ] : "6 
Sa oe Soe ee 
sin Z5 sin Z 5 sin 25 g sin a 
3.2707100 «7726223 -5110397 = 0414955 
»7726223 10.3315155 4,8919810 = 1753744 
9110397 4, 8919810 3.7896887 = 0940528 


Solution of the above normal equations gives in minutes of are 


Aa, = 0'.009 
Z ' 

Ag, = 0",013 
digs OG 

A¥q = 0'.007 


These quantities show no improvement is justified. 


5. Application of the Gravity Direction Cosines to a Gravity Referenced 
Datum 


There are a number of possibilities. The simplest notion is to allow 
the intersection of the 1%, plane with the Xo Yo plane to be the Yo axis 
and the intersection of the Xe Y. plane with the Xo Yo plane be the Xo axis. 


This concept is illustrated in Figure 6. 
Given 


COS A, COs Bas COS Yq 
tan Fel ae 
cos BG 
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Figure 6 SIMPLEST CONCEPT OF GRAVITY REFERENCE PLANE 
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whence 


cos «o = cos 6 


Gy G 
cos PGs aus 86 
cos MGs = 0 
cos "Gy = sin Tel cos Yq 
cos Ba = cos Tal cos Yq 
Y 
cos Ye, = sin Yo 


Such a datum orientation is not likely to be entirely satisfactory 
inasmuch as the arbitrary datum selected is probably inherent to the object 


such as the X axis defining the axis of the object. In such a case, it may 


be desirable to let Xo or Yo be the preferred axis projected to the Xa Yo 


plane from Za This selection is illustrated in Figure 7. Since [Z, Xo = 90° 


cos ag = sin AG 
. xX 


from the low of cosines 


0 = cos a 


C cos Ba + sin a 


sin Ba cos s 


G RY 


0 = cos Ba cos Y, + sin Be sin Yg 008 Sy, 


0 = cos Ye cossa + sin'y, sin a, cos s 


G G ZX 
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Figure 7 PROJECTION OF ARBITRARY AXIS Xqg TO Xg IN PLANE 2Z¢Xq 
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whence 
iiica =. cos AG cos Ba 
ae sin sin 8 
tamara © 
Kgercs = 2 608 8. COS Yq 
ate sin sin 
G G 
Resta = OS Yq COs o 
a sin sin 
G 
Therefore 
cos Ba = -sin Ba aa el oe 8G 
xX sin o, sin 8g 
G G 
cos Ya, = -sin y ee. I Geese SG 


The direction of Yo is deduced from 


cos a 
G. 


¥ a 


cos el cos el 


¥: xX 


cos a = 
Gy sin 
cos Ge = cos Y, (sin Ot) = COS Og (-cos y 
2 
= cos ¥¢ (sin Ge 4 cos*a 
sin el 


ite! 


so Gee 
sin Ye sin el 
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= cos By (-cos Ye cos % ) = (-cos 8 


G 


cos a 


G 


sin a 


cos a 


G 


sin al 


cos © 


sin a 


COS Yq 


G 


cos the + cos By cos Bit cos Y¢ cos te* 0 
Y 


+ cos g cos 8 + cos y cos y = 0 
G G G G 
x Y x 


G ) cos Yo * 0 
G 


) 
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COS Yq 
sin Oe 
cos y = cos a, (-cos B Seay - cos 6. sin a 
ey : eras. aie 
G 
cos Be 
sin el 


Summarizing for the projection of either Xs or Yo. 


a 
x Projected qe Projected 
cos o = sin o cos oO eg 
G G G : 
Xx Xx sin G 
cos Bo = =cos £6 Pee cos Bg =Q 
Xx sin o xX 
reer cos Y¥_ = COs 4 
cos y, = = cos y G xX sin Bo 
xX sin o 
=e . cos B 
es a 0 cos ce cos oa G 
Y sin ¢ 
cos 8 a) pee Ye cos 8B = + sin 8 
G : G G 
Y sin o ¥ 
G 
cos y, =- cos +, ~°8 8 
cos a, SOF Ba: Gy sin 8. 
Ye : sin Ba 
Y sin o 


G 


In the most general sanse, any preferred line 6 may be projected to 
the Xo Yo plane. Assume it is desired to project any selected line to the 
Xo Ya plane and let the projected line be the Xo or Yq axis. The direction 


cosine of Zo and o are known whence 
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cos cos ~ cos g. cos. 
cos a ve Yo BS Yo 


x sin G 
fo) 


cos Yo: cos a. .7. COS .Y, cos a 


cos B = G 


X sin G 
fe) 


cos a,.cos B. —- cos.a. cos... 
¢ G By ° we 


COS Yq 
X sin Gy 

where 

cos Go = cos O cos a, + cos Be cos BY + cosyy, cos y, 
Then 

cos "Gy = cos g cos Ye, - cos Pe, COS Yq, 

cos Pay = COS Yq COs cm ~ cos Yo, cos el 

cos Ye, ; cos we cos a, - cos *e, cos Ba 


The projection of any line is illustrated in Figure 8. If it is 
desired to let Yq be the perpendicular to the Zo plane, only the subscripts 


are changed. 


In any case, whatever the decision the matrices and coordinates are 


rotated to a gravity reference: 


given 
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Figure 8 PROJECTION OF ANY LINE 0 TO THE Xg Yg PLANE 
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cos 8 cos Ba cos. 8B 


cos cos Y, cos ¥ 
G. G. G 
xX XY 


and an orientation matrix 


cos a cos O cos a 

x y Zz 
cos cos cos 

BS BY BE 
cos. cos cos 

¥- tg Ve 


The gravity reference matrix is reduced as follows: 


= + + 
cos O, Cos a, COS dq cos BY cos BG cos Y, COS Yq 
& X x X 
cos a = cosirq cos a + cos B_ cos 8 + cos cos 
y. y Gy ¥ G. Vy ue 
& 4 
cos G | = COS 4 COS OG + cos 8. cos Ba + cos Y, ©O8 Yq 
: Xx X x 
= + 
cos B. cos a, cos Jel + cos Be cos em cos Vg cos for 
g x 


cos 8 =cos a cosa, +t cos ® cos B + cos y_ cos y 
Vp ¥ Gy ¥ Gy y Gy 
cos Be = cos a, cos Ge + cos 6. cos ene + cos Vo S98 Ye. 
g 
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= + 
cos Y, cos a, cos oe cos B. cos Ba + cos Y,, cos Ye 
& 
cos = cos a cos a, + cos B cos 8. + cos cos 
ae y G y G Yy fe 
cos ‘a, = cos a, cos al + cos B. cos Ba + cos Y, cos Ye 


More important perhaps is the rotation of the space coordinates to the 


gravity referenced datum 


i] 
ba 


X cos a + Y cos 8 + Z cos y 
a Gy a Gy a Gy 


i] 
ky 


X cos a + Y_ cos g + 2 cos y 
a Gy a Gy a Gy 


x cos Og + YT, cos BG + Zs cos Yo =Z 
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